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Abstract. The Multidimensional Assignment Problem (MAP or s-AP 
in the case of s dimensions) is an extension of the well-known assignment 
problem. The most studied case of MAP is 3-AP, though the problems 
with larger values of s have also a number of applications. In this paper 
we propose a memetic algorithm for MAP that is a combination of a 
genetic algorithm with a local search procedure. The main contribution 
of the paper is an idea of dynamically adjusted generation size, that yields 
an outstanding flexibility of the algorithm to perform well for both small 
and large fixed running times. The results of computational experiments 
for several instance families show that the proposed algorithm produces 
solutions of very high quality in a reasonable time and outperforms the 
state-of-the art 3-AP memetic algorithm. 

Keywords: Multidimensional Assignment Problem, Metaheuristic, Ge- 
netic Algorithm, Memetic Algorithm. 



1 Introduction 

The Multidimensional Assignment Problem (MAP or s-AP in the case of s di- 
mensions) is an extension of the well-known Assignment Problem (AP, linear AP) 
which is exactly two dimensional case of MAP. MAP has a host of applications, 
see, e.g., [1, 3,4] for 'classic' applications and [2, 12, 13] for recent applications in 
solving systems of polynomial equations and centralized multisensor multitarget 
tracking. In fact, several applications described in [2, 3, 12, 13] naturally require 
the use of s-AP for values of s larger than 3. 

For a fixed s > 2, s-AP is stated as follows. Let Xi = X2 = . . . = Xg = 
{1, 2, . . . , 71} and let X = XiX X2 x . . . x Xg- For a vector e G X, the component 
Cj denotes its jth coordinate, i.e., ej G Xj. Each vector e G X is assigned a non- 
negative weight w{e). A collection A oi t < n vectors e^, e^, . . . , e* is a (feasible) 
partial assignment if 7^ e'j holds for each i ^ k and j e {1,2,..., s}. The 
weight of a partial assignment A is w{A) = X)i=i 'w{e''). An assignment (or full 
assignment) is a partial assignment with n vectors. The objective of s-AP is to 
find an assignment of minimum weight. 

While AP can be solved in a polynomial time [10], s-AP for every s > 3 
is NP-hard [5]. MAP is a very hard problem in the following sense. The weight 



matrix of MAP contains values, there exist n!*~^ possible assignments and the 
fastest known algorithm to find the optimal one takes 0(n!''~^n^) time. Indeed, 
without loss of generality set = i for i = 1, 2, . . . , n and find the optimal values 
of ej, by solving a linear AP for every possible combination of the values e*- for 

1 = 1, 2, . . . , n and j = 2, 3, . . . , s — 1 such that ej- ^ Cj for i ^ k. Compare it 
with, e.g., the Travelling Salesman Problem which has only weights, (n — 1)! 
possible tours and that can be solved in 0(n^2") time. 

The rest of the paper is organized as follows. The general scheme and all the 
details of the proposed memetic algorithm including the local search procedure 
are described in Section 2. The test bed for our computational evaluation is 
discussed in Section 3. The experiments with our and some other metaheuristics 
are described and the results are discussed in Section 4. The main outcomes of 
the presented research are reported in Section 5. 

2 The Algorithm 

A memetic algorithm is a combination of genetic algorithm with local search. A 
typical scheme of a memetic algorithm is as follows. 

1. Produce the first generation, i.e., a set of solutions. 

2. Apply the local search procedure to every solution in the first generation. 

3. Repeat the following while a termination criterion is not met: 

(a) Produce a set of new solutions by applying so-called genetic operators 
to solutions from the previous generation. 

(b) Improve every solution in this set with the local search procedure. 

(c) Select several best solutions from this set to the next generation. 

While the general scheme of the algorithm is quite common for all memetic 
algorithms, the set of genetic operators and the way they are applied can vary 
significantly. In our algorithm, we use the following procedure to obtain the next 
generation: 

gf'+i = selection{{gl} U mutation{g'' \ {g{},Pm, IJ-m) U C), 

where g'^ is the fcth generation and g'- is the jth. assignment of the fcth generation; 
g\ is the best assignment in the fcth generation. Constants pm = 0.5 and Hm = 
0.1 define the probability and the strength of mutation operator respectively. 
The function selection simply returns m,+i best distinct assignments among 
the given ones, where rrife is the size of the fcth generation (if the number of 
distinct assignments in the given set is less than rrii+i, selection returns all the 
distinct assignments and updates the value of mi+i accordingly). The set of 
assignments C (crossover part) is calculated as follows: 

(p-mi+i-mi)/2 

C = Local S ear ch{cr ossover {g\^, g\j)) 



where u,v S {1,2, . . . ,mi} are chosen randomly for every crossover run and 
p = 3 defines how many times more assignments should be produced for the 
selection operator. The mutation function for a set of solutions is defined as 
follows: 



where r G [0,1] is chosen randomly every time. The functions crossover {x,y), 
perturb{x, /u) and Local S ear ch{x) are discussed later. 

2.1 Coding 

In genetic algorithms, coding is a way to represent a solution as a sequence of 
atom values such as boolean values or numbers; genetic operators are applied to 
such sequences. Good coding should meet the following requirements: 

~ Coding should be a bijcction, i.e., a decoding procedure such as decode {code (x)) = 
X for any feasible solution x should exist. 

— Evaluation of tlic^ quality of a coded solution should be fast. 

— Every fragment of the coded solution should refer to just a part of the whole 
solution, so a small change in the coded sequence should not change the 
whole solution. 

— It should be relatively easy to produce algorithms for random modification 
of a solution (mutation) and for combination of two solutions (crossover) 
which produce feasible solutions. 

Huang and Lim [7] use a local search procedure that, having first two dimen- 
sions of an assignment, determines the third dimension (note that the algorithm 
from [7] is designed only for 3-AP). Since the first dimension can always be 
fixed with no loss of generality, one needs to store only the second dimension 
of an assignment. Unfortunately, this coding requires a specific local search and 
is robust for 3-AP only. We use a different coding: a vector of an assignment is 
considered as an atom in our algorithm and, thus, a coded assignment is just a 
list of its vectors. The vectors are always stored in the first coordinate ascending 
order, e.g., an assignment 



2.2 Termination Condition 

Usually, a termination condition in a memetic algorithm tries to predict the 
point after which any further effort is useless or, at least, not efficient. A typical 




{(2,4,1),(4,3,4),(3,1,3),(1,2,2)} 



would be represented as 



(1,2, 2), (2, 4,1), (3, 1,3), (4, 3, 4). 



approach is to count the number of subsequent generations which did not im- 
prove the best result and to stop the algorithm when this number reaches some 
predefined value. 

Wc use a different approach. To be able to compare different algorithms cor- 
rectly and to satisfy real world requirements, we bound our algorithm within 
some fixed running time. Apart from the mentioned advantages of this termina- 
tion condition, it is worth to note that it gives flexibility to produce either fast 
or high quality solutions depending on one's needs. 

2.3 Generation Size 

The most natural way to fit the running time of a memetic algorithm into the 
given bound is to produce generations of some fixed size until the time is elapsed. 
However, it is clear that one memetic algorithm cannot work efficiently in both 
cases if there are just a few generations and if there are hundreds of generations. 
Thus, instead, we fix the number of generations and vary the generation size. 

Our computational experiments show that, with a fixed running time, the 
most appropriate number I of generations for our algorithm is always around 50; 
this number does not depend on the local search procedure or the given time. 
Since the running time of the local search procedure can vary significantly (e.g., 
last generations usually contain better solutions than the first ones and, thus, 
are processed faster) and also to make our algorithm easily portable, we decided 
to adjust the generation size dynamically according to the remained time such 
that the total number of generations would always be close to /. 

In particular, the size of the next generation is calculated as follows: 



where T is the given time, t is the elapsed time, A is the time spent to produce 
the previous generation, I is the prescribed number of generations and k = 1.25 
is a constant that limits the generation size change. The 'otherwise' branch 
normally should not be used but in case of quick finish of the Ith generation 
it increases the size of generation (recall that the goal is to keep the number 
of generations close to /). Note that the values are real numbers, and the 
actual size rrii of the ith generation is defined as a rounded down value of . In 
addition, the value p ■ rrij+i — rrii, i.e., the number of assignments produced by 
crossover, should always be even because every crossover produces exactly two 
assignments. If this condition is not met, the value of mj+i is increased by one: 



The size of the first generation is obtained in a different way (see Section 2.4). 





2.4 First Generation 



As it was shown in [6] (and we also confirmed it by experimentation with our 
memetic algorithm and the construction heuristics from [8]), it is good to start 
any MAP local search or mctaheuristic from a Greedy construction heuristic. 
Thus, we start from running the Greedy algorithm (we use the same implemen- 
tation as in [6]) and then perturb it using our perturb procedure to obtain every 
item of the first generation: 

gj = Local S ear ch{perturb{gr eedy, fif)), 

where greedy is an assignment obtained by the Greedy heuristic and /x/ = 
0.2 is the perturbation strength coefficient. Since perturb performs a random 
modification, it guarantees some diversity in the first generation. 

The algorithm produces assignments for the first generation until T/I time 
elapses or at least 4 assignments (recall that T is the time given for the whole 
memetic algorithm and / is the prescribed number of generations). 

2.5 Crossover 

A typical crossover operator combines two solutions, parents, to produce two 
new solutions, children. Crossover is the main genetic operator, i.e., it is the 
source of a genetic algorithm power. Due to the selection operator, it is assumed 
that good fragments of solutions are spread wider than others and that is why, 
if both parents have some similar fragments, these fragments arc; probably good 
and should be copied without any change to the children solutions. Other parts 
of the solution can be randomly mixed and modified though they should not be 
totally destroyed. 

The widest used crossover is the one-point crossover that produces two chil- 
dren x' and y' from two parents x and y as follows: x- = x, and = yi for 
every i = 1,2, ... ,k and x'^ = yi and y[ = Xi for every i = k + 1, k + 2, . . . ,n, 
where k e {1,2, . . . ,n — 1} is chosen randomly. One can see that if Xi = yi for 
some i, then the corresponding values in the children sequences will be the same: 

However, the one-point and some other standard crossovers do not preserve 
feasibility of MAP assignments since not every sequence of vectors can be de- 
coded into a feasible assignment. We propose a special crossover operator. Let x 
and y be the parent assignments and x' and y' be the child assignments. First, 
we retrieve equal vectors in the parent assignments and initialize both children 
with this set of vectors: 

x' = y' = X fly . 

Let k = \x n y\, i.e., the number of equal vectors in the parent assignments, 
p = x\x' and r = y\y', where p and r are ordered sets. Let tt and u be random 
permutations of size n — k. For every j = 1, 2, . . . , n — fc the crossover sets cither 
x' = x' U p^(j) and y' = y' U r^Q) or x' = x' U r^^j) and y' = y' U p^^y The 
particular option is chosen randomly with the probability of the first one 80%. 



Since this procedure can yield infeasible assignments, it requires additional 
correction of the child solutions. For this purpose, the following is performed for 
every dimension d = 1,2, ... ,s for every child assignment c. For every i such 
that 3j < i : cjj^ — set = r where r £ {1,2, . . . ,n} \ {c\, c^, . . . , c^} is chosen 
randomly. In the end of the correction procedure, the vectors in the assignment 
should be sorted in the ascending order of the first coordinates (see Section 2.1). 

In other words, our crossover copies all equal vectors from the parent assign- 
ments to the child ones, then copies the rest of the vectors randomly choosing 
every time a pair of vectors, one from the first parent and one from the second 
one, and then adding them either to the first and to the second child respectively 
or, with probability 20%, to the second and to the first child respectively. Since 
the obtained child assignments can be infeasible, the crossover corrects each of 
them; for every dimension of every child it replaces all duplicate coordinates with 
randomly chosen correct ones, i.e., with the coordinates which are not currently 
used for that dimension. 

2.6 Perturbation Algorithm 

The perturbation procedure perturb{x, ji) is intended to modify randomly an 
assignment x with the given strength /x. In our memetic algorithm, perturbation 
is used to produce the first generation and to mutate assignments from the 
previous generation when producing the next generation. 

Our perturbation procedure pertor6(x, ^) performs [n/Lt/2] random swaps. 
Each swap selects two vectors and some dimension randomly and then swaps 
the corresponding coordinates: swap and y^, where u,v G {1, 2, . . . , n} and 
d G {1, 2, . . . , s} are chosen randomly; repeat the procedure [n/x/2] times. For 
example, if = 1, our perturbation procedure modifies up to n vectors in the 
given assignment. 

2.7 Local Search Procedure 

An extensive study of a number of local search heuristics for MAP can be found 
in [6]; the paper includes both fast and slow algorithms. It also shows that a 
combination of two heuristics can yield a heuristic superior to the original ones. 

The following heuristics were considered as possible candidates for the local 
search procedure for our memetic algorithm (we provide only a brief description 
of every heuristic here; full descriptions can be found in [6]): 

— 2-opt is a simple heuristic that selects the best of all possible interchanges 
for every pair of vectors in the assignment. 2-opt is known as a very fast but 
poor quality heuristic. 

— DV is a natural extension of a local search proposed and used in [7]. On every 
iteration, DV fixes all but one dimensions yielding a 2-AP instance, solves 
it (recall that 2-AP can be solved in a polynomial time) and reorders the 
unfixed dimension according to the obtained solution of 2-AP. 



— MDV is a modification of DV proposed in [6]. While DV divides the dimension 
set into {1,2,..., . . . s} and {i} for every i, MDV tries every possible 
division, fixes dimensions from every of two sets together yielding a 2-AP 
and then proceeds as DV. For 3-AP, DV and MDV are equivalent while for 
,s-AP with ,s > 3 MDV is more powerful. 

— DV2 is a combination of 2-opt and DV; it applies sequentially 2-opt and DV 
to the given assignment. DV2 repeats the whole procedure until one of the 
heuristics appears to be incapable to improve the assignment. 

— MDV2 is like DV2 but instead of DV it apphos MDV. 

~ MDV3 is like MDV2 but instead of 2-opt it applies 3-opt. 3-opt, similarly to 
2-opt, tries every possible interchange for every triple of vectors in the as- 
signment. While 3-opt is a very slow and not very successful heuristic, MDV3 
is relatively fast and outperforms all other considered local search heuristics 
with respect to solution quality. 

— MDVV is like MDV2 but instead of 2-opt is applies v-opt. v-opt is an exten- 
sion for s-AP with arbitrary s of the Variable Depth Interchange heuristic 
proposed in [1] for 3-AP. Like 2-opt, v-opt considers interchanges of vector 
pairs, however the objective and the enumeration order in v-opt are totally 
different. 

Results for 3-opt and v-opt as a local search for our memetic algorithm are 
not given in this paper since they did not show any promising results in our 

experiments. 

Table 1 shows comparison of the results of our memetic algorithm based 
on different local search procedures. One can see that the fast heuristics MDV, 
DV2 and MDV2 [6] perform better than others in almost every experiment, and 
MDV2 shows the best average result among them. Similar results were obtained 
for the Multichain metaheuristic proposed in [6] (the results reported in [6] are 
slightly different but it can be explained by the difference in the test beds) . The 
situation stays approximately the same in the experiments with different given 
times. Thus, MDV2 was selected as a local search procedure for our memetic 
algorithm. 

3 Test Bed 

In this section we discuss instance families used for experimental evaluation of 
our memetic algorithm. We use pscudo real-world instances: Composite Cyclic 
(CC), Composite Clique (CQ), Squares Root (SR) and perturbed versions of each 
of them. The CC and CQ families were proposed before in the literature (see [3] 
and references there) while the others arc new. 

To explain these instance families, it is better to represent MAP as a graph 
problem. Let G = {Vi U V2 U . . . U Vg, £■) be an s-partite graph where |Vi| = n 
for every i = 1, 2, . . . , s. A weight is assigned to every clique in G, i.e., to every 
induced subgraph {{vi, V2, Vs}, { V1V2, V1V3, viVg, V2V3, V2V4, V2Vs, 
. . . , Vs-iVg}), where vi £ Vi, V2 £ V2, . ■ . , Vg & Vg-, a. clique corresponds to a 



vector in the original interpretation. The objective is to find n disjoint cliques 
such that their total weight is minimized. 

The idea of all CC, CQ and SR is that in the real life it is usually possible to 
define some relation only between two objects but not between s objects. Thus, 
the graph G should bo weighted and the weight of a clique C = {Vc, Ec) should 
be a function of the weights Eq of the edges in the clique. We immediately get 
two natural instance families, CC and CQ: 



The SR instance family, as well as CC, is based on the weights in a cycle: 



In this case, the objective is not only to minimize the considered weights but 
also to keep all the weights not too large. 

The perturbed instances are the same instances but the weight assigned to 
every clique is modified randomly: w(c) = w{c) + r for every clique c, where 
r G {0, 1, ... , 19} is chosen randomly. 

Our test bed includes instances of 3-AP, 4-AP, 5-AP and 6-AP, six types 
of instances per each value of s: CC, perturbed CC, CQ, perturbed CQ, SR and 
perturbed SR (for s = 3, CQ and perturbed CQ were excluded from the test 
bed since CC and CQ for 3-AP are equivalent; however, the results for the 3-AP 
CC (perturbed CC) instances are included in the CQ (perturbed CQ) averages). 
The sizes n of the instances are as follows: n = 40 for 3-AP, n = 30 for 4-AP, 
n = 18 for 5-AP and n = 12 for 6-AP. Those sizes were obtained empirically such 
that the difficulty of all instances in the test bed is approximately the same (in 
other words, a local search takes approximately equal times for every considered 
instance). We produced 10 different instances for every experiment, i.e., for every 
combination of s, n and instance family; the results reported in the tables in 
Section 4 in every row are averages among 10 runs. To generate the instances, we 
initialized every edge of the graph G with a random number from {1.2,..., 100} 
using standard Miscrosoft .NET random generator [11] which is based on the 
Donald E. Knuth's subtractive random number generator algorithm [9]. As a 
seed of the random number sequences for all the instance types we use the 
following number: seed = s + n + i, where i € {1,2,..., 10} is an index of the 
instance of this type and size. 

4 Experiment Results 

Three metaheuristics were considered in our experiments: 
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— A Multichain metaheuristic (MC) proposed in [6] which significantly outper- 
forms all local search heuristics known from the literature with respect to 
solution quality. 

— An extended version of the memetic algorithm by Huang and Lim [7] (HL). 

— Our memetic algorithm (GK). 

For the purpose of comparison, the Huang and Lim's heuristic was extended 
for s > 3 as follows: 

— The coded assignment contains not only the second dimension but it stores 
sequentially all the dimensions except the first and the last ones, i.e., an 
assignment {e^, e^, . . . , e-'} is represented as el, e^, . . . , 62 , e|, e§, . . . , 63 , 
. . . , el_i, el_i, . . . , e"_i {e\ = 1 for each i and e\ can be chosen in an 
optimal way solving an AP). 

— The local search heuristic, that was initially designed for 3-AP, is extended 
to DV as it is described in [6]. 

— The crossover, proposed in [7], is applied separately to every dimension (ex- 
cept the first and the last ones) since it is designed for one dimension only 
(recall that the memetic algorithm from [7] stores only the second dimension 
of an assignment, see Section 2.1). 

Our computational experience show that the solution quality of our imple- 
mentation of the Huang and Lim's heuristic is similar to the results reported 
in [7] and the running time is slightly larger because of the extension for s > 3. 

All the heuristics are implemented in Visual C-|— I- and evaluated on a plat- 
form based on an AMD Athlon 64 X2 3.0 GHz processor. 

The results are reported in Tables 2 and 3. The following running times are 
considered in this tables: 0.3 s, 1 s, 3 s, 10 s and 30 s. The first given time, 
0.3 s, corresponds approximately to the nmning time of HL if it applies its 
original termination condition. Every entry of these tables contains an average 
solution error for 10 instances of some fixed type and size but of different seed 
values (see Section 3 for details) . The value of the solution error is calculated as 
(y — Ubest)/'f^best ' 100%, where v is the obtained solution value and Vhest is the 
best known solution value^ (the best solution values can be found in Table 1). 
The instance name consists of throe parts: the number of dimensions s, the type 
of the instance and the size n of the instance. In addition, letter 'p' is added to 
the instance name to indicate that it is a perturbed instance. 

The average values for different instance families and numbers of dimensions 
are provided at the bottom of each table. The CQ (perturbed CQ) averages 
include CC (perturbed CC) results for 3-AP since CC and CQ instance families 
are equivalent when s — 3 and, thus, CQ and perturbed CQ instances for 3-AP 
are excluded from the test bed. The best among MC, HL and GK results are 
underlined in every row for every particular given time. 

One can see that GK clearly outperforms MC and HL: GK performs better than 
others in most experiments both with small given times (0.3 and 1 seconds) and, 

^ The best known solutions were obtained during our experiments with different 
heuristics. 



especially, with larger given times (3, 10 and 30 seconds) for which it is the best 
in every experiment with just a few exceptions. GK shows the best average results 
for every instance family, for every number of dimensions s and for every given 
time. The latter fact proves its ability to perform well for very different given 
times. The MC heuristic mostly outperforms HL for small given times but shows 
approximately the same or slightly worse results for larger given times that was 
expected because memetic algorithms are normally intended to be high quality 
but slow heuristics. 

It is worth to note that we experimented with different values of the GK 
algorithm parameters such as I, fx f, fim etc. and concluded that small variations 
of these values do not significantly influence the algorithm performance. 

5 Conclusion 

A memetic algorithm for MAP is proposed in this paper. The main feature of 
the algorithm is a dynamically adjusted generation size. It was shown that there 
exists the most appropriate number of generations which does not depend on the 
given instance or solution quality requirements. This idea allowed us to design 
a flexible algorithm that is able to perform well for very different given times. 
Extensive experiments show that the proposed algorithm clearly outperforms all 
high quality MAP algorithms known from the literature given the same running 
time. 
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Table 1. Comparison of the memetic algorithm based on different local search. The 
given time is 1 sec. 
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0.94 


0.17 


0.15 


0.10 


0.05 


3.30 


1.72 


5cql8p 3617.0 


1.16 


0.24 


0.60 


0.14 


0.24 


3.58 


0.94 


5srl8 504.9 


4.08 


3.37 


0.85 


0.87 


0.34 


8.69 


1.86 


5srl8p 610.5 


4.73 


4.46 


2.62 


2.24 


0.84 


12.33 


3.03 


6ccl2 1023.8 


0.24 


0.79 


0.00 


0.17 


0.12 


3.87 


0.43 


6ccl2p 1108.3 


0.17 


0.59 


0.10 


0.13 


0.03 


3.46 


0.22 


6cql2 4505.6 


0.28 


0.16 


0.26 


0.00 


0.30 


3.51 


1.17 


6cql2p 4608.9 


0.16 


0.07 


0.31 


0.00 


0.22 


2.69 


0.51 


6srl2 502.9 


0.50 


0.68 


0.02 


0.08 


0.04 


9.33 


0.46 


6srl2p 567.8 


0.46 


1.02 


0.35 


0.18 


0.05 


5.30 


1.95 


All avg. 


4.96 


2.45 


1.35 


1.07 


0.66 


6.58 


2.47 


CC avg. 


5.12 


2.45 


1.10 


0.89 


0.59 


6.15 


2.14 


CC p. avg. 


4.63 


2.68 


1.40 


1.16 


0.44 


6.02 


2.61 


CQ avg. 


3.60 


0.88 


0.78 


0.46 


0.57 


4.50 


1.90 


CQ p. avg. 


3.17 


0.71 


0.86 


0.45 


0.46 


4.45 


1.42 


SR avg. 


7.65 


3.58 


1.61 


1.12 


0.72 


8.65 


2.58 


SR p. avg. 


6.86 


3.64 


2.10 


2.07 


1.08 


8.50 


3.37 


3-AP avg. 


8.15 


0.77 


0.95 


0.47 


0.65 


4.62 


0.98 


4-AP avg. 


9.97 


5.77 


3.18 


2.71 


1.53 


9.27 


5.78 


5-AP avg. 


2.50 


2.17 


0.96 


0.79 


0.32 


7.07 


1.82 


6-AP avg. 


0.30 


0.55 


0.17 


0.09 


0.13 


4.69 


0.79 



Table 2. Metaheuristics comparison. 



Solution error, % 





l).:5 sec. 






1 S(>C. 






•'i sec. 




Inst. 


HL 


MC 


GK 


HL 


MC 


GK 


HL 


MC 


GK 


OCC4U 


2.48 


2.18 


i.04 


u.yi 


L.oZ 


n /I 


U.oU 




U.JO 


oCC4Up 


2.25 


1.32 


1 1 r 


U.4e 


U.y / 


n /1 
U.4o 


U.44 


u. < y 


r\ A 
U. J4 


osr4U 


z .oo 


9 ^9 


Z.Uo 


u.y / 


i.o < 


U.4 1 


r\ AQ 




n on 


3sr40p 


3.62 


2.84 


2.19 


0.88 


2.02 


1.09 


0.58 


1.56 


0.23 


4cc30 


13.38 


2.24 


3.06 


3.82 


1.30 


1.41 


1.92 


0.86 


0.50 


4cc30p 


13.16 


4.43 


4.35 


5.47 


3.50 


1.06 


2.03 


2.60 


0.62 


4cq30 


6.39 


3.02 


2.29 


1.44 


2.58 


1.41 


0.90 


2.30 


0.54 


4cq30p 


5.91 


2.73 


2.33 


1.58 


1.90 


0.91 


1.10 


1.36 


0.73 


4sr30 


16.69 


4.72 


4.18 


8.25 


3.06 


2.04 


3.29 


1.74 


2.20 


4sr30p 


18.95 


6.57 


4.71 


7.71 


4.16 


2.35 


2.94 


3.18 


0.93 


5ccl8 


0.34 


1 /II 

i.4i 


0.80 


1.72 


0.89 


0.29 


1.40 


0.65 


0.06 




7.13 


2.01 


9 IT^ 


2.22 


1 fin 


n 1 8 


1 Qfi 


1 fin 


n 98 


5cql8 


0.68 


1.68 


0.35 


0.25 


1.57 


0.05 


0.10 


1.57 


0.02 


5cql8p 


0.61 


1.50 


0.63 


0.29 


1.14 


0.24 


0.23 


0.67 


0.09 


5srl8 


8.66 


1.56 


1.19 


2.97 


0.81 


0.34 


2.40 


0.81 


0.36 


5srl8p 


10.19 


3.72 


2.69 


3.62 


3.60 


0.84 


2.47 


3.60 


0.29 


6ccl2 


1.05 


1.16 


0.34 


0.75 


1.05 


0.12 


0.70 


0.49 


0.00 


6ccl2p 


2.17 


0.22 


0.29 


1.10 


0.22 


0.03 


0.76 


0.22 


0.00 


6cql2 


0.21 


1.21 


0.83 


0.21 


0.99 


0.30 


0.20 


0.93 


0.03 


6cql2p 


0.10 


0.89 


0.28 


0.10 


0.86 


0.22 


0.07 


0.86 


0.02 


6srl2 


3.50 


0.52 


0.36 


1.53 


0.22 


0.04 


1.31 


0.18 


0.00 


6srl2p 


5.64 


1.44 


1.02 


1.90 


0.74 


0.05 


1.90 


0.67 


0.00 


All avg. 


6.00 


2.27 


1.76 


2.19 


1.67 


0.66 


1.26 


1.35 


0.36 


CC avg. 


5.81 


1.75 


1.44 


1.80 


1.21 


0.59 


1.13 


0.84 


0.23 


CC p. avg. 


6.18 


1.99 


1.96 


2.32 


1.57 


0.44 


1.30 


1.30 


0.31 


CQ avg. 


2.44 


2.02 


1.25 


0.70 


1.69 


0.57 


0.42 


1.54 


0.23 


CQ p. avg. 


2.22 


1.61 


1.10 


0.61 


1.22 


0.46 


0.46 


0.92 


0.30 


SR avg. 


7.93 


2.33 


1.94 


3.43 


1.49 


0.72 


1.86 


1.12 


0.71 


SR p. avg. 


9.60 


3.64 


2.65 


3.53 


2.63 


1.08 


1.97 


2.26 


0.37 


3-AP avg. 


2.81 


2.21 


1.73 


0.81 


1.62 


0.65 


0.49 


1.38 


0.30 


4-AP avg. 


12.41 


3.95 


3.49 


4.71 


2.75 


1.53 


2.03 


2.01 


0.92 


5-AP avg. 


5.60 


1.98 


1.28 


1.85 


1.60 


0.32 


1.43 


1.48 


0.18 


6-AP avg. 


2.11 


0.91 


0.52 


0.93 


0.68 


0.13 


0.82 


0.56 


0.01 



Table 3. Metaheuristics comparison. 



Solution error, % 







10 sec. 






.'iO s(>c. 




Inst. 


HL 


MC 


GK 


HL 


MC 


GK 






l.UO 


n 99 
u.zz 


u.oy 


u. / o 


n 09 
u.uz 


occ4Up 


U.zo 


U.I / 


U.U 1 


n 0/1 
U.z4 


U.I 1 


n no 


osr4U 




1 /I o 

1.4z 


U.Uj 




i.Uo 


U.UU 


3sr40p 


0.41 


1.36 


0.10 


0.31 


1.21 


0.02 


4cc30 


1.55 


0.53 


0.23 


1.43 


0.39 


0.10 


4cc30p 


1.57 


1.73 


0.47 


1.43 


1.62 


0.19 


4cq30 


0.85 


1.70 


0.52 


0.74 


1.47 


0.11 


4cq30p 


0.88 


1.05 


0.09 


0.88 


0.89 


0.06 


4sr30 


3.04 


0.84 


0.99 


2.74 


0.58 


0.11 


4sr30p 


2.53 


1.63 


0.21 


1.73 


1.39 


0.37 


5ccl8 


1.00 


0.65 


0.00 


0.85 


0.65 


0.00 




1 Sfi 
± .ou 




n on 


± .uu 


1 97 


n nn 

U.UU 


5cql8 


0.04 


1.57 


0.04 


0.04 


1.57 


0.02 


5cql8p 


0.18 


0.67 


0.07 


0.17 


0.67 


0.00 


5srl8 


1.96 


0.81 


0.16 


1.96 


0.81 


0.00 


5srl8p 


2.26 


3.60 


0.13 


1.75 


3.36 


0.00 


6ccl2 


0.70 


0.40 


0.00 


0.70 


0.32 


0.00 


6ccl2p 


0.62 


0.22 


0.00 


0.62 


0.22 


0.00 


6cql2 


0.20 


0.93 


0.00 


0.20 


0.93 


0.00 


6cql2p 


0.04 


0.86 


0.00 


0.04 


0.86 


0.00 


6srl2 


1.11 


0.04 


0.00 


0.95 


0.04 


0.00 


6srl2p 


1.78 


0.53 


0.00 


1.66 


0.30 


0.00 


All avg. 


1.07 


1.09 


0.15 


0.94 


0.96 


0.05 


CC avg. 


0.92 


0.66 


0.11 


0.84 


0.52 


0.03 


CC p. avg. 


1.08 


1.08 


0.13 


0.99 


0.97 


0.07 


CQ avg. 


0.38 


1.31 


0.19 


0.34 


1.17 


0.04 


CQ p. avg. 


0.34 


0.84 


0.06 


0.33 


0.80 


0.04 


SR avg. 


1.60 


0.78 


0.30 


1.48 


0.63 


0.03 


SR p. avg. 


1.74 


1.78 


0.11 


1.36 


1.56 


0.10 


3-AP avg. 


0.34 


1.15 


0.10 


0.30 


0.95 


0.03 


4-AP avg. 


1.74 


1.25 


0.42 


1.49 


1.06 


0.16 


5-AP avg. 


1.22 


1.48 


0.07 


1.07 


1.39 


0.00 


6-AP avg. 


0.74 


0.50 


0.00 


0.70 


0.44 


0.00 



